The spontaneous emergence of ordered phases in crumpled sheets 
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X-ray tomography is performed to acquire 3D images of crumpled aluminum foils. We develop an 
algorithm to trace out the labyrinthian paths in the three perpendicular cross sections of the data 
matrices. The tangent- tangent correlation function along each path is found to decay exponentially 
with an effective persistence length that shortens as the crumpled ball becomes more compact. In 
the mean time, we observed ordered domains near the crust, similar to the lamellae phase mixed 
by the amorphous portion in lyotropic liquid crystals. The size and density of these domains grow 
with further compaction, and their orientation favors either perpendicular or parallel to the radial 
direction. Ordering is also identified near the core with an arbitrary orientation, exemplary of the 
spontaneous symmetry breaking. 
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Crumpling is capable of producing a highly rigid struc- 
ture with a record minimum of material. Every child 
knows how to make a baseball out of a crumpled news- 
paper. Even with the assistance of a dutiful parent, its 
interior shall remain roughly 80% vacant pQ. In addition 
to this application, equally fascinating and puzzling phe- 
nomena related to crumpling abound in a wide range of 
length scales, for instance, how DNA is packed in the 
tiny capsule of viruses [2 or the practical challenge for 
the auto industry to come up with a better design to 
protect the safety of drivers during car accidents [3]. 

The enormous resistance of a crumpled ball can be 
attributed to the geometrical constraint and the self- 
avoidance. The former refers to the inevitable develop- 
ment of D-cones due to the unstretchability of a thin 
sheet [U [5]. However, beyond the geometric and mechan- 
ical properties of a single and two-cone interaction [6], 
collective behavior of the microstructures such like ridges 
and vertices remains unexplored. Previous simulations [7] 
have identified the phantom and self- avoiding sheets as 
belonging to different universality classes because the 
universal exponent a of their force-size relation is dif- 
ferernt. To clarify the effect of self- avoidance, x-ray to- 
mography becomes highly desired because it enables us 
to study the evolution of the internal structure system- 
atically and perform calculations with its data [8] [9]. 

It is known [4 that the macroscopic properties are 
shaped more by the collective behavior rather than in- 
dividual ridges and vertices. This is because these singu- 
larities are linked by not just the tensile force and bend- 
ing rigidity, but also correlated by the strong hard core 
interact ion [7| [10]. When a flat sheet is subject to a gen- 
tle force, the first deformation due to the buckling is of 
a conical shape [TT]. As the compaction progresses, the 
single cone deformation is followed by a large number 
of ridges and vertices [1, while the facets they encircle 
also begin to align. In this Letter, we present the first 
systematic analysis of the local and global structural or- 
dering inside the crumpled aluminum ball. An algorithm 
is developed to trace out the curves in x-ray tomography 



and reconstruct the cross-section view without destroy- 
ing the sample. Through the buckling and ordering, we 
study how they accumulate and affect the final configu- 
ration in this highly non-Markov process. 

Nine sheets of aluminum foil with different diameter 
(Ro[mm] = 3, 4, 5, 6, 6.5, 7, 8, 9, 10) are randomly folded 
by hand first and then squeezed by the flat tip tweezer at 
different directions [12] into balls of the same final radius 
R =1.5 mm. To determine whether they still belong 
to the thin sheet regime, we calculate their Foppl-von 
Karman number[4], (Rl/h 2 )[12 x (1 — v 2 )] where v = 0.35 
is the Poisson ratio and h = 16//m is the thickness of the 
foil. Ranging from 2.4 x 10 4 to 2.6 x 10 5 , they turn out 
to be of the same order as in previous work[TJ [7J fTU| [TTj . 

A special version of microtomography is employed, 
based on the high intensity x-ray from synchrotron[13 . 
It provides a standard resolution between 1-2 /am which 
shows clear reconstructed images for our analysis. The 
experiment is performed at the 01 A beamline of National 
Synchrotron Radiation Research Center in Taiwan. The 
beamline provides unmonochromatic x-rays whose en- 
ergy distribution is 8-15 keV. Image acquisition time per 
projection is about 10ms, which is captured by a CCD 
with 2X optical lens focused on a CdW04 single crys- 
tal scintillator. The resulted reconstruction consists of a 
data matrix of 1200x1200x1200 pixels of size 3 fim. 

To study the packing configuration quantitatively, it 
is necessary to vectorize the data points. This segmen- 
tation method is complex and case-dependent [14 . The 
most challenging part in tracing the crumpled surface 
is to distinguish two contact planes. According to previ- 
ous work [15] which concluded that cross sections through 
different angles share the same statistics, we are assured 
that each of these cross sections is representative of the 
bulk configuration. Therefore, we can focus on the devel- 
opment of 2D tracing algorithm and present its result as 
a precursor for a full 3D construction from these images. 

After all the images have been reconstructed into 3D 
data matrices, we resample three perpendicular cross sec- 
tions which are X-Y, Y-Z, Z-X planes of the sample. 
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Briefly, the procedures of tracing algorithm are mainly 
divided into three parts: random seeding, identifying so- 
lutions, and labeling traced points. We start with a circle 
around the seed point with an radius of 24 pixels, and 
then find all crossing points between its perimeter and 
the paths. Since many paths can be close to each other, 
multiple solutions frequently occur. When this happens, 
we select the point that is joined to the seed point from 
all candidate solutions. The traced points are marked 
immediately to avoid double tracing with the forest-fire 
simulation. After all the points have been vectorized, we 
perform a high order Beizer fitting to resample the seg- 
mented points to reduce irregularities and increase the 
sampling rate, see Fig{l] A more detailed description 
of the tracing algorithm and segmented images can be 
found in the online supplemental document [14J. 




FIG. 1: Panel (a) shows a slice of raw images reconstructed 
from 1000 projections for R/Ro — 0.167. The glisterns and 
rings are experimental artifacts [5] that can be reduced by the 
fill tracing method. As a contrast, the segmented configura- 
tion is shown on panel (b) by linking the traced points. 
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FIG. 2: Tangent- tangent correlations C(As) are plotted for 
three different compactions R/Ro. Although the data fluctu- 
ates in the tail, which is likely an artifact of too few samplings 
of long paths, an exponential decay and a persistence length 
can be extracted. Thick lines illustrate the experimental data, 
and thin lines are the fitting curves. The inset highlights the 
decrease of persistence length as the ball shrinks. 

We start by calculating the tangent-tangent correla- 
tion function which is a basic statistical property related 



to the buckling [16] of sheets. This function is defined 
as C(s,s') = ((u(s) - (u(s))) • (u(s') - (u(s')))), where 
u(s) = dH/ds is the tangent vector at arc length s in the 
curvilinear coordinate and R(s) is the position vector. It 
can be simplified to C(\s — s'\) = C(As) if the system 
exhibits translational invariance. This assumption was 
checked quantitatively to hold except at the end points 
where the fluctuations become large. 

The function C(As) is found to decay exponentially in 
Fig j2] with an effective persistence length l p . This decay 
form can be derived from the random packing of facets. 
Borrowing the concept of tube model for polymers [17], 
one can think of the sheet as moving inside two walls 
which model the confinement due to the hard core inter- 
action from its neighboring portions of sheet. When we 
cut perpendicularly through the walls, the cross section 
will reveal a wiggling path with a static configuration 
similar to that of a polymer in the tube model. However, 
it should be noted that their dynamics are different be- 
cause the path is, afterall, a projection of a 2D sheet. The 
movement of each segment needs to coordinate with the 
rest of the sheet, unlike the reptation model in which the 
polymer is confined in all sidewise directions by a static 
tube. The spacer width apparently decreases as we in- 
crease the crumpling force. Dividing this width by the 
segment length a of the path gives the maximum angle 
between neighboring segments. Roughly, we can imag- 
ine the configuration as being mapped out by a random 
walk with a fixed stride a but only two choices of angular 
deviations, ±0. Then the probability of finding the rela- 
tive angle between the nth and zeroth segments equaling 
mO obeys the Gaussian distribution for a random walk 
after n = (s f — s)/a steps: P(m, n) ~ exp(— m 2 /2n). The 
function C(As) = (cos(ra#)) can be explicitly evaluated: 




where l p = 2a /0 2 . The inset of Fig J2] shows that l p short- 
ens as R/Ro decreases, which implies the crumpled struc- 
ture becomes more disordered along the curvilinear di- 
rection. This observation also requires that the segment 
length a not only depends on the bending rigidity of the 
material, but can be cut short by the compact packing. 

The concept of an effective persistence length is simi- 
lar to that in polymers [17] and the de Gennes coherence 
length in membranes [1 8 . However, different from the 
latter two cases, our crumpled ball is so macroscopic [19] 
that the thermal temperature becomes irrelevant. In- 
stead, it is the noise from the random folding that allows 
the sheet to appear rumpled. This zero-temperature ran- 
domizing effect also exists in the granular systems [20]. 

After calculating the order along the curvilinear coor- 
dinate, we turn to the cartesian plane. Note that locally 
ordered structures can be identified in Fig|l] The paths 
in the plane, which cut through the facets in the 3-D 
sample, show a tendency to align and form lamella-like 
phase. We can separate the ordered portion from the 
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FIG. 3: We trace out trajectories and define different domains 
marked by quadrangles for R/Ro = 0.15. Different color or 
gray level signalizes different orientations. The domain forma- 
tion caused by the hard core interaction and high density of 
paths is similar to that of lyotropic liquid crystals [2T], 
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FIG. 4: Four properties of the lamella-like phase are calcu- 
lated as a function of R/Ro. They include the number of 
lamellae in panel (a), ratio of total area and mass covered by 
the lamellae in (b), size of each lamella in (c), and mass or to- 
tal length encompassed in each lamella in (d). Data in (a) and 
(d) can be fit by 1.57 (R/Ro)- 1 - 05 and 10 5 exp(-7.04R/#o). 



disordered one to define crystalline and amorphous re- 
gions. This classification is aided by the vectorized data. 
In Fig(3j all paths are denoted by different colors and 
brightnesses in grayscale to indicate different orientation 
to the (1,0) direction. The neighboring facets with the 
same orientation are marked with the same color, and 
can be easily recognized as a lamella. To extract the do- 
main boundaries, the Laplacian of the color brightness 
was first calculated with the local maxima signalizing 



the location of the divisions. Although the distribution 
of lamellae can in principle be derived by this method, 
manual identifications are still required when the color 
gradient is too noisy. Therefore, we use the quadrangles 
to label all lamellae with the reference of Laplacian field 
to obtain the data with a coherent format. 

Using the data averaged over three perpendicular cross 
sections for each sample, we calculate four essential prop- 
erties of the lamellae: their number, ratio of total area 
and mass they covered in the cross section, their size, and 
the mass encompassed in each lamella. We checked that 
the product of data in panels (a) and (c) equaled the area 
ratio in (b) times the cross section area, ttR 2 . Figure [I] 
shows that they all increase as R/Ro decreases. Panel (a) 
indicates the number of lamellae is inversely proportional 
to R/Ro. Two features are worth noting in (b): Firstly, 
the reason why the two lines are not plainly proportional 
is that the total mass grows as we fix R and increase Ro 
to achieve lower R/Ro. In contrast, total area ttR 2 is 
unchanged. Furthermore, the alignment allows for more 
efficient packing inside the domain which explains why 
the mass ratio is larger. Secondly, both ratios never ex- 
ceed 0.5. This is similar to supercooled liquid where the 
extent of ordering is hindered from being complete [22J. 
One may wonder how a structure with so much amor- 
phous region can be so hard. A possible explanation is 
that these ordered domains near the boundary interlock 
and act like a hard crust. Each domain consists of many 
aligned layers which greatly enhance the bending rigidity 
and make them more resistant to buckling. 




Compaction R/Ro 

FIG. 5: The bulk order S is shown to grow when the sample 
becomes more compact. The solid line is a guide to the eyes 
to emphasize the trend that the ordering increases as the ball 
becomes more compact. The i?o(mm)=3, 4, and 5 data are 
not inlcuded because they no longer exhibit a shell region and 
also their error bars are too large due to insufficient samplings. 
The inset shows the wall-rod correlation functions G(r) for 
R/Ro =0.25 and 0.19. 

In the previous work [23 that demonstrates a sponta- 
neous patterning in vibrated rods, a wall-rod correlation 
function G(r) = (cos(2rj)) was defined to extract the size 
of the core or bulk area, where () denotes averaging over 
(f and r] is the angle between the tangent vectors at point 
(r, ip) and the boundary in the ip + ir/2 direction. We 
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repeat the same definition to study the packing configu- 
ration near the core and plot the results in the inset of 
FigjEJ Since the external force acts from outside, it is 
natural that the correlation with the boundary decreases 
as we enter the core. The division between the bulk and 
shell regions is marked by the first vanishing of G(r). 
Again, following the notation of [23], an order parameter 
£=(cos(20)) is defined for the bulk where 9 is the an- 
gle between the tangent vector and lamella direction of 
the bulk. We are surprised to find a spontaneous bulk 
ordering in Fig(5] where correlations with the boundary 
layer have considerably weakened. As R/Ro decreases, 
the bulk order increases. The ordered phase induced by 
high concentration is similar to that happened in [23] and 
lyotropic liquid crystals. According to Onsager[24], al- 
though parallel arrangements of anisotropic objects lead 
to a decrease in orientation entropy, there is a gain in 
positional entropy. Thus, a positional order is expected 
to become entropically favorable at sufficient rod concen- 
trations. 

Given that the deformations in aluminum are plastic 
and irreversible, crumpling can be viewed as a series of 
quenching process since not all configurations are acces- 
sible, nor equally probable. The noise introduced by the 
random folding still enables the sheet to slightly adjust its 
configuration to seek a local potential minimum [25 . This 
random process plays the role of vibration in the sponta- 
neous patterning of vibrated rods [23] and slow shearing 



in granular systems [20j |26]. Since this large amount of 
facets are not only correlated by the sheet but also inter- 
act strongly via the bending potential, the geometrical 
constraint forbids our crumpled sheet from reaching the 
true ground state [27]. Consequently, the bulk order pa- 
rameter of our compact sample is still much lower than 
that of vibrated rods [23 . 

In conclusion, we performed x-ray tomography to 
quantitatively study the inner structure of aluminum foils 
at different compactions. All paths in the three perpen- 
dicular cross sections are traced out and vectorized before 
calculating the statistical properties of their packing con- 
figuration. The tangent-tangent correlation of the path 
reveals an effective persistence length that decays with 
the compaction. A second length scale associated with 
the size of domains that emerges near the crust and mim- 
ics the lamella phase in lyotropic liquid crystals. Number 
of these domains and their area and mass grow monoton- 
ically with the compaction, and their orientation favors 
either perpendicular or parallel to the radial direction. 
We also identified an ordered domain near the core of 
the crumpled ball with an arbitrary orientation, exem- 
plary of the spontaneous symmetry breaking. 
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